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Observing a physical quantity without disturbing it is a key capability for the control of individual 
quantum systems. Such quantum non-demolition (QND) measurements were first introduced in the 
1970s in the context of gravitational wave detection to measure weak forces on test masses by 
high precision monitoring of their harmonic motion. At present, QND techniques have become an 
indispensable tool in quantum science for preparing, manipulating, and detecting quantum states 
of light, atoms, and other quantum systems. Here we experimentally perform QND measurements 
of the position of a mechanical oscillator by using short pulses of light, i.e. with a pulse duration 
much shorter than a period of mechanical oscillation. This technique enables us to perform both 
state preparation and full state tomography of the mechanical motional state. We have reconstructed 
mechanical states with a position uncertainty reduced to 19 pm, limited by the quantum fluctuations 
of the optical pulse, and we have performed 'cooling by measurement' to reduce the mechanical mode 
temperature from an initial 1100 K to 16 K. Future improvements to this technique may allow for 
quantum squeezing of mechanical motion, even from room temperature, and reconstruction of non- 
classical states exhibiting negative regions in their phase-space quasi-probability distribution. 



Experiments are now beginning to investigate non- 
classical motion of massive mechanical devices [T]-[5]. 
This opens up new perspectives for quantum-physics en- 
hanced applications and for tests of the foundations of 
physics. A versatile approach to manipulate mechanical 
states of motion is provided by the interaction with light, 
typically confined to microwave or optical cavities. Such 
cavity-optomechanics experiments j6HTQ] have thus far 
largely concentrated on high sensitivity continuous mon- 
itoring of the mechanical position pTHT6] . Because of 
the back-action imparted by the probe onto the measured 
object, the precision of such a measurement is fundamen- 
tally constrained by the standard quantum limit [17]rT8]. 
and therefore only allows for classical phase-space re- 
construction [m [191 HO] • In order to observe quantum 
mechanical features that are smaller than the mechani- 
cal zero-point motion, back-action-evading measurement 
techniques are required. Following the early insights of 
Braginsky, beating the standard quantum limit 'can he 
achieved only in one way: design the probe so it ''sees'' 
only the measured observable' ^7\. Such back-action 
evading techniques were first realized for the detection 
of optical quadratures [2TJ [22] and have since been used 
for, e.g. precision measurement of atomic ensemble spin 
[23H27] and QND microwave photon counting [28 . In 
optomechanics, to perform a back-action evading mea- 
surement of the mechanical position, a time-dependent 
measurement scheme is required. One prominent exam- 
ple is the so-called 'two-tone approach' [29l [30], which 
uses a probe with an intensity that oscillates at twice the 
mechanical frequency. The field probes the mechanics 
periodically and the back-action imparted to the mechan- 
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ical motion by the optical probe does not affect the mea- 
surement of the mechanical amplitude. This is closely 
analogous to a stroboscopic measurement of the mechan- 
ical motion [29 . Using the two-tone approach with a mi- 
crowave probe field, a back-action evading measurement 
of a single quadrature of a nanomechanical resonator has 
recently been demonstrated [31] . 



Here we take a different tack by using single optical 
pulses instead of a periodic interaction. Our experiment 
employs optical pulses that have a duration much shorter 
than a mechanical period of motion. This enables a back- 
action-evading measurement of the mechanical position 
because the interaction leaves the position unchanged, 
perturbing only the mechanical momentum, and was first 
suggested by Braginsky [32]. We implement the pulsed 
protocol proposed in Ref. [33] where one or two pulses 
are used to prepare a motional state 'by measurement' 
and then a subsequent pulse is used for state tomogra- 
phy. The precision of this measurement process is ulti- 
mately limited by the optical quantum noise, which we 
have reached in our experiment. We therefore require a 
quantized description of the optical field, however, at this 
stage, the mechanical motional state may be described 
classically. Utilizing a pulsed protocol has the signifi- 
cant advantage that the experiment can be performed on 
a timescale faster than decoherence or rethermalization 
and thus has considerable tolerance to the surrounding 
thermal bath [33]. Also, employing our pulsed approach, 
mechanical dynamics rather than the steady-state can 
be conveniently probed and non-equilibrium mechanical 
behavior can be characterized. Also note that pulsed 
quantum optomechanics operates fully in the so-called 
'non-resolved sideband regime', in which the cavity de- 
cay rate is much larger than the mechanical frequency. 
Indeed, all results reported here were obtained without 
the use of an optical cavity. 



Experimental protocol 

Our experimental setup is shown schematically in 
Fig. [ija). Optical pulses are injected into a Mach- 
Zehnder interferometer that has a micro-mechanical os- 
cillating mirror in one of the two interferometer paths. 
The pulses are first divided by a beam-splitter that forms 
one intense beam that acts as a local oscillator (LO) and 
one weak beam that we will henceforth refer to as the sig- 
nal. The signal is focussed onto and reflects from a micro- 
mechanical oscillator (Fig. [ijb)). During the reflection 
of the short optical pulse, changes to the position of the 
mechanical oscillator are negligible. The coherent optical 
pulse gains a phase shift in proportion to the mechani- 
cal position, which is accurately described by a phase 
quadrature displacement as the mechanical position fluc- 
tuations are small. Concurrently, the radiation-pressure 
force of the reflection imparts momentum to the mechan- 
ical resonator. This momentum can be decomposed into 
a classical component due to the mean photon number 
and a component dependent upon the photon number 
fluctuations. Quantitatively, this optomechanical inter- 
action is described by the input-output relations 
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Here, the subscripts label the light (L) and mechanics 
(M); X and P are the dimensionless amplitude (posi- 
tion) and phase (momentum) quadratures for the light 
(mechanics); x = ^^XQ\fN /\ quantifies the quadrature 
information exchanged between the light and the me- 
chanics and determines the strength of the mechanical 
position measurement, and Vt = SttxqN/X is the classi- 
cal momentum transfer to the mechanical oscillator. (A/", 
mean photon number per pulse; A, optical wavelength; 
xq = {h/2mQf^ujMY^'^ mechanical ground state exten- 
sion; meff, mechanical effective mass; uom^ mechanical 
angular frequency.) After the optomechanical reflection, 
the signal then overlaps and interferes with the LO pulse 
on a 50/50 beam-splitter where the (mean) phase be- 
tween the LO and signal beams is set to be 7r/2. The in- 
tensities of both beam-splitter outputs are measured by 
photodiodes and the photocurrents are subtracted to im- 
plement homodyne detection of the optical phase quadra- 
ture. A typical difference current time trace is shown in 
Fig. [TJc) where the measurement outcome Pl is the time 
integral over the pulse duration of the difference current. 

After the pulsed measurement the mechanical state of 
motion is changed as our knowledge of the mechanical po- 
sition has increased. For an initial thermal state of the 
mechanical resonator with a large thermal occupation, 
i.e. x^(l + 2n) > 1, the means and variances of the me- 
chanical quadratures, upon obtaining the measurement 
outcome Pl are [33 : 
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FIG. 1: (a) Schematic of the experimental setup used to per- 
form state tomography and state preparation of the motional 
state of a mechanical resonator. In addition to the optical 
pulses a weak continuous field is used to stabilize the inter- 
ferometer phase using the homodyne output passed through 
a low-pass filter with cutoff frequency below the mechanical 
frequency, (b) Colourized optical micrograph of the high- 
reflectivity micro-mechanical cantilever fabricated for this ex- 
periment. The head of the cantilever, where the signal beam 
is focussed, is 100 /xm in diameter, (c) Example time trace 
of the homodyne output for a pair of 4 /is pulses. (For clar- 
ity, the pulse rising and falling edges are not shown.) The 
measurement outcome Pl is the time integral of the homo- 
dyne output (indicated by the shaded region). Time resolved 
optical quantum noise is visible during the pulse. 



where n ^ kBT/hwu is the mean occupation of the 
mechanical mode when in thermal equilibrium with the 
environment at temperature T > hwM/^B- Notably 
the measurement reduces the mechanical position vari- 
ance from n to 1/(2%^) , which does not depend on 
the initial occupation. The resultant state of mechani- 
cal motion, following such a measurement, is no longer 
in thermal equilibrium with the surrounding environ- 
ment and has a reduced effective thermal occupation 
neff = - 1/2 ^ (n/(2x2))i/2. Moreover, 

a subsequent pulse performed after one quarter of me- 
chanical harmonic evolution can measure the mechanical 
momentum at the time of the first pulse to further reduce 
the effective occupation. This 'cooling-by-measurement' 
method for entropy reduction is rapid and has consider- 
able tolerance to the bath temperature [33 . With future 
experimental improvements, this scheme allows for the 
generation of high purity and quantum squeezed states 
of mechanical motion 'by measurement'. Due to the re- 
silience to mechanical thermal noise, this scheme may 
provide a more feasible route to quantum squeezing than 
achieving a quantum squeezed state via parametric mod- 
ulation [11] |19] , which can be combined with continuous 
measurement and feedback [34]. 
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In our experiment one or two pulses are used to pre- 
pare a mechanical state at a known time. Follow- 
ing state preparation a read-out pulse is made after 
time 9/ujm of mechanical harmonic evolution to sam- 
ple the mechanical probability distribution of the 6- 
rotated quadrature, i.e. a marginal. Repeating this 
process many times and obtaining the marginals for 
a large number of mechanical phase-space angles is 
sufficient to uniquely determine the mechanical quan- 
tum state of motion [35 . Quantum state tomogra- 
phy by measurement of the marginals was ffist realized 
with optical fields using homodyne interferometry [5B] 
and has now become an indispensable tool in the field 
of quantum optics [37] being applied to other physical 



systems such as molecular vibration [38^, spin ensem- 
bles [39 , and microwave fields [40^. Here we implement 
such mechanical state tomography [33] by utilizing the 
pulsed measurement outcome probability distribution 
Pr(Pi) = /dXM7r-i/2exp[-(Pi - xXM)']Pr(XM,0) 
that contains the mechanical marginals Pr(XM,^) = 
{Xm \ P^Mi^) I^m), where p^^{0) is the mechanical input 
state to be reconstructed after time 9/uom of harmonic 
evolution. In this experiment we prepare and reconstruct 
mechanical motional states with features that are not 
smaller than and hence, unless otherwise noted, we 
use the optical measurement outcome distribution as an 
approximation for the mechanical distribution using the 
scaled outcome Pl/x !33] . 



Mechanical motional state reconstruction and state 
preparation via measurement 

The mechanical resonator used for this experiment is 
a micro-mirror cantilever constructed from an epitaxial 
AlxGai_xAs crystalline multilayer, see Fig. 1(b). The use 
of such a monocrystalline material structure allows for a 
significant reduction of the mechanical damping of the 
resonator when compared with dielectric refiectors [4T] 
and simultaneously provides high optical refiectivity. The 
crystalline material used here is nominally identical in 
composition and individual layer thickness to the struc- 
tures in Ref. [42 and is designed for maximum refiectiv- 
ity at our optical wavelength of 1064 nm. The multilayer 
mirror comprises 40.5 layer pairs in order to minimize 
transmission losses. The cantilever was etched from a 
6.88 /im thick multilayer and is 1.45 mm in length with a 
cantilever arm 5 /im in width with a circular head 100 /im 
in diameter where the optical signal beam is focussed [43] . 
Note that the resonator is etched directly from the mul- 
tilayer mirror material and is therefore equally refiec- 
tive at all points along the structure with an (inten- 
sity) reflectivity of 99.982 %. This cantilever has a fun- 
damental out-of-plane vibrational mode with frequency 
co'M/27r = 984.3 Hz, effective mass meff = 260 ng [44] . 
and a mechanical quality of Q = 3.1 x 10^ in vacuum 
(10~^ mbar) and at room temperature measured via me- 
chanical ringdown. 

Our optical setup (Fig. 1(a)) was constructed from 
optical-fiber-based components that provided good phase 
stability and excellent spatial mode matching. Indeed, 
when the optical powers in the two arms of the interfer- 
ometer are balanced we observed an interference visibility 
exceeding 99.9 %. We use a continuous laser source and 
generate optical pulses of duration 1 /is (excluding the 
pulse edges) with a fiber-based intensity modulator. The 
mean photon number in a signal pulse was up to 10^ and 
in order to provide a homodyne signal well above the elec- 
tronic noise we use a large LO to signal ratio with up to 
10^^ photons per LO pulse. (These photon numbers were 



determined via optical power measurement during con- 
tinuous wave operation.) The signal pulses are directed 
onto the cantilever head using an anti-refiection coated 
fiber focuser and are then retro-refiected. To calibrate 
the proportionality between the measurement outcomes 
and the mechanical position we refiect the signal beam 
from a rigid mirror adjacent to the mechanical resonator 
and scan the mirror position using a calibrated piezoelec- 
tric actuator recording both the piezo scan positions and 
measurement outcomes [44] . 

After a pulse measurement is performed to sample a 
mechanical marginal, the mechanical state is reinitial- 
ized by first allowing it to return to equilibrium with 
the environment and then the mechanical state is re- 
prepared. This process is repeated many times to accu- 
mulate sufficient data to characterize the statistical prop- 
erties of the mechanical motion. The marginal distribu- 
tions were then obtained by constructing a histogram 
from the many measurement outcomes recorded for each 
mechanical phase-space angle 0. As the states studied 
here are symmetric about the Xm and Pm axes we mea- 
sure a set of many marginals with angles between = 
and = 7r/2 to fully characterize the state of motion. 
The phase-space probability distribution W{Xm^Pm) is 
then obtained by using the inverse Radon transformation 
on the set of marginals. 

The measurement results we obtained for motional 
state preparation and reconstruction are summarized in 
Fig. |2j In Fig. |2|a) a reconstruction of an initial thermal 
state that is driven by white noise across the mechanical 
resonance up to a mode temperature of 1100 K is shown. 
This temperature was obtained using the equipartition 
theorem ksT^f^ = mef^ujj^cr'^^ where the mechanical po- 
sition variance was determined from the calibrated 
measurement outcome distribution after subtracting the 
optical noise contribution. A single pulsed measurement 
made on this initial thermal state generates a motional 
state that has a reduced position uncertainty (Fig.[2|b)). 
The momentum distribution of this state is negligibly 
changed as the back-action to the mechanical momen- 
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FIG. 2: Motional states of a mechanical resonator experimentally prepared and fully reconstructed using pulsed quantum 
measurement. The uppermost row shows the pulsed protocols used; pulses filled pink are used for measurement based state 
preparation and pulses filled red are used for mechanical state tomography. On the two rows below, a subset of the mea- 
sured marginals and the reconstructed phase-space distributions are shown, respectively. The phase-space distributions were 
reconstructed using 9 marginal angles up to ^ • 7r/180 = 90° (with a larger number of bins used than that shown for the 
marginals). To help provide a sense of scale, the marginal distributions are plotted with mechanical displacement in meters 
and the phase-space distributions are plotted in units of the mechanical zero-point motion, (a) Tomography and reconstruction 
of an initial mechanical thermal state driven by white noise up to a mode temperature of 1100 K. The dashed circle has a 
radius equal to 2a of the thermal distribution, (b) A single pulsed measurement reduces the mechanical position variance, 
but leaves the momentum distribution unchanged, (c) 'Cooling by measurement' performed with two pulses separated by one 
quarter of a mechanical period rapidly reduces the mechanical state's entropy. The effective temperature of the mechanical 
state reconstructed here has been reduced to 16 K. (d) State reconstruction of a non-Gaussian mechanical state of motion 
generated by resonant sinusoidal drive, (e) The (one standard deviation) width of the position distribution observed for states 
(a-c) with phase-space angle 0. (The left vertical axis is in units of nanometers and the right axis uses units of mechanical 
quantum noise.) The thermal state (red points) shows a position width approximately twice that when at room temperature 
(dashed line). State (b) has a reduced position width for small phase-space angles (purple points). The position width of state 
(c) is reduced for all phase-space angles (blue points). The solid lines are theoretical fits obtained using Eq. ^ generalized 
for all as well the two-pulse-preparation case (see Ref. [33]). (f) Conditional mechanical width with pulse strength measured 
using preparation and read-out pulses separated by 5° of mechanical evolution. The dashed line is a theoretical fit with a model 
using the two units of optical quantum noise and the finite mechanical evolution. The solid line is the inferred conditional 
mechanical width immediately after the preparation pulse. The vertical line indicates the pulse strength used for states (a-c). 



turn made by the reflection of the optical pulse is much with a random but known mean due to the random mea- 
smaller than the large initial momentum uncertainty. surement outcome, see Eq. ([2|. By making the transfor- 
Each pulsed measurement generates a mechanical state 
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(r) and (pi) indicate read-out and preparation, respec- 
tively, this random mean is subtracted and the distribu- 
tion of the mechanical state can be characterized. We 
would like to emphasize here that no 'post-selection' is 
performed and all measurement outcomes are used in this 
process. The plots for Figs. ^di,h) were generated from 
the same data set where the statistics of the prepara- 
tion pulse alone characterizes the unconditional initial 
thermal state and the read-out pulse characterizes the 
conditional mechanical state. 

In Fig. [2|c) the reconstruction of a mechanical state 
of motion prepared via two pulsed measurements sepa- 
rated by one quarter of a mechanical period is shown. 
For this pulse sequence the read-out pulse outcome is 
transformed using P^'^^ P^j^^ - P^j^^^ cos 9 + P^^'^ sin l9, 
where is the angle of mechanical evolution made be- 
tween the second preparation pulse and the read-out 
pulse. The width of the mechanical phase-space distribu- 
tion has been significantly reduced in both the position 
and momentum quadratures compared to the initial ther- 
mal state (Fig. [2|a)) and hence the effective mode tem- 
perature has decreased. Ideally, for this mechanical state, 
the width of the mechanical marginals should be constant 
for all ^, however, in our experiment the phase corre- 
lation between the pulses reduces with increasing pulse 
separation as low frequency noise, due to imperfect phase 
locking, enters the signal. This results in a broadening of 
the conditional mechanical marginals as 9 increases. The 
effective temperature Teff = mQf^uj\^G^~^^ /ks 
observed for this state is 16 K, which depends on the 
product of the standard deviations of the position and 
momentum quadratures. Were the pulses to remain cor- 
related to within the quantum noise, the effective tem- 
perature that could be reached for this measurement 
strength, taking the effects of mechanical rethermaliza- 
tion into account is 4.4 K [33 . We would like to highlight 
here that rethermalization contributes to only less than 
1 % of this value. To summarize the observed effects 
of single and two-pulse mechanical state preparation dis- 
cussed above. Fig. [2|e) provides a plot of the measured 
mechanical widths with 9 for the initial thermal state 
and the two mechanical conditional states. In this plot 
the mechanical widths were determined from the cali- 
brated pulse outcome distributions after subtracting the 
optical noise contributions that were measured indepen- 
dently. The data for both of the mechanical conditional 
states was taken with the same signal pulse powers and 
for each phase space angle 300 pulses were recorded to 
construct the histograms. 

As an example of a non-Gaussian state of motion we 
have reconstructed a driven thermal state (Fig.[2|d)) that 
was generated by applying a sinusoidal drive on reso- 
nance with the mechanical eigenfrequency. Note that 
the two peaks in the mechanical marginals are narrower 
than the broad thermal state in Fig. 2|a) as this state was 
prepared at room temperature without the white noise 
drive. Even though this state of motion and the ther- 



mal state are rotationally invariant in phase-space many 
marginals are measured for their reconstruction. On the 
other hand, the conditional mechanical states of motion 
are not rotationally invariant in phase-space as the time 
of the preparation pulse (s) sets the time for 6> = 0. Note 
that this pulse-based tomography scheme does not mea- 
sure the angle ^ = as the read-out pulse is temporally 
separated from the preparation pulse (s). The lack of this 
marginal angle causes the rippling near Xm = in the 
reconstructed phase-space distributions. By employing 
shorter pulses and measuring the marginals at smaller 
angles this rippling can be easily reduced. 

To demonstrate the scaling of our measurement 
strength in Fig. |2|f) the conditional mechanical width 
observed by a read-out pulse made after 5° of mechani- 
cal free evolution is plotted with increasing pulse ampli- 
tude. For this pulse separation the two pulses are well 
correlated and the width of the conditional mechanical 
state is limited by the optical quantum noise in the mea- 
surement (see the methods section for more details). As 
the signal pulse strength is increased the standard devi- 
ation of the conditional mechanical position distribution 
decreases with A/""^/^, which is a result of the optical 
number-phase uncertainty relation. The dashed line in 
the plot is a theoretical prediction including the two units 
of optical shot noise, one each for the preparation and 
read-out pulses, and the small contribution from the me- 
chanical evolution between the two pulses. The relative 
amplitudes for the data points were measured precisely 
and scaled by a free fitting parameter into units of square- 
root photon number, where the photon number per pulse 
obtained is consistent with measurements of the opti- 
cal power made during continuous wave operation. For 
the largest optical pulse strength used the statistics of 
the read-out pulse demonstrate a conditional mechanical 
width (after the preparation pulse) of = 19 pm corre- 
sponding to a measurement strength of x = 2.1 x 10~^. 



Discussion 

The techniques developed in this work provide the 
ability to experimentally perform quantum optomechan- 
ics in the time domain. This offers significant poten- 
tial for optomechanics-based quantum information and 
quantum metrology applications by providing the frame- 
work for quantum state preparation of a mechanical res- 
onator via quantum measurement [45]. One may then 
also envision combining such measurement based state 
preparation with feedback to implement full quantum 
control [46]. One exciting example of mechanical dy- 
namics that can be probed by pulsed optomechanics has 
been recently theoretically discussed by Buchmann et 
al. [47 , where pulsed measurements, as now realized in 
this work, are considered for the observation of quan- 
tum tunneling of a mechanical oscillator in a double- well 
potential. Another example for quantum state prepa- 
ration is that, even though the optomechanical inter ac- 
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tion used here is linear with the mechanical position, by 
exploiting the optical non-linearity, measurements 
with a strength significantly larger than that attainable 
with dispersive optomechanics can be performed [48 . An 
measurement can be used to conditionally prepare 
highly non-Guassian mechanical superposition states and 
experimentally characterizing the decoherence of such 
states is important to determine the feasibility of using 
mechanical elements for coherent quantum applications 
and can also be used to empirically test collapse mod- 
els [49-52]. The pulsed measurements performed here 
may also be utilized for a QND-measurement-based light- 
mechanics quantum interface [53 . Furthermore, a se- 
quence of four pulsed optomechanical interactions can 
be used to generate non-classical mechanical states of 
motion via an optomechanical geometric phase [54 and 
can even be used to experimentally explore potential 
quantum-gravitational phenomena [55 . 

For the preparation of quantum squeezed states of me- 
chanical motion the measurement strength needs to be 
increased to x > 1- An effective route to meet this re- 
quirement would be to employ an optical cavity to en- 
hance the optomechanical interaction. Using the experi- 
mental parameters achieved in this work, a cavity finesse 
of 10^ is sufficient. As such a cavity simultaneously re- 
quires a high finesse, as well as a large bandwidth to 
accommodate a short optical pulse, this is best achieved 
with an optomechanical microcavity [33]. Such improve- 
ments to the measurement sensitivity will not only en- 
able Wigner reconstruction with significant negativity 
but, owing to this pulsed protocol's resilience against 
mechanical thermal noise, may also allow the generation 
of non-classical mechanical states in the regime of room 
temperature quantum optomechanics. 



Methods 

To verify that the measurement scheme used here is 
optical quantum noise limited we measured the phase 
quadrature conditional variance of a pair of optical pulses 
with increasing total photon number, i.e. the sum of the 
signal and LO photons per pulse, while keeping the sig- 
nal to LO ratio fixed, see Fig. [3] As with our calibration 
procedure, the signal beam is focussed onto a rigid mirror 
adjacent to the mechanical oscillator to prevent coupling 
to the mechanical motion. The pulse separation used for 
this measurement was 14.1 /is, which would correspond 
to 5° of mechanical free evolution, and is the same as 
that used for the dataset shown in Fig. [2|f). With this 
pulse separation the conditioning is essentially the sec- 
ond pulse outcome minus the first pulse outcome. The 
quantum noise components of these two temporal modes 
are uncorrelated, however, the lower frequency classical 
noise components vary slowly between the two pulses and 
are thus suppressed by the conditioning. Quantum me- 
chanics predicts a linear dependence for the variance with 
total photon number, whereas, were classical phase noise 



to be the dominant contribution, a quadratic dependence 
with the total photon number per pulse would be ob- 
served [56 . During this measurement we were limited 
to a total photon number of 10^^ as the phase lock per- 
formance dramatically reduced beyond this point. Were 
we able to measure beyond this optical power the clas- 
sical phase noise would eventually become the dominant 
noise and the conditional mechanical variance that can 
be achieved would saturate. 

The data points for Fig.[3]were obtained from Gaussian 
fits to histograms of the conditional outcomes. The error 
bars indicate a one standard deviation uncertainty as de- 
termined from the fit. The observed conditional variance 
shows a linear dependence with the total photon num- 
ber with a 'goodness of fit' parameter = 0.97, taking 
the error bars into account. This demonstrates that, up 
to a total photon number of order 10^^, the conditional 
variance is quantum noise limited. 

Also included in Fig.jsjis the measured electronic noise, 
i.e. the conditional variance observed using no light. This 
contribution is 19.5 dB smaller than the observed optical 
quantum noise at the data point with the highest optical 
intensity (A^tot = 9-5 x 10^). 
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A^TOT [xlO^] 

FIG. 3: The measured optical phase quadrature conditional 
variance (2) (i) plotted with the total photon number (sig- 

nal plus local oscillator) per pulse Atot- The linear depen- 
dence observed demonstrates that the measurement scheme is 
optical quantum noise limited up to Atot — 10^°. Were clas- 
sical phase noise to be the dominant noise source a quadratic 
dependence would be observed. 
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Supplementary Information - Experimental Pulsed 
Quantum Optomechanics 

Effective mass measurement 

The optically probed effective mass of a mechanical vi- 
brational mode depends upon (i) the geometry and mate- 
rial properties of the mechanical structure and (ii) the in- 
tensity profile of the incident optical beam. The mass as- 
sociated with the mechanical displacement mode shape, 
i.e. the modal mass, is in general less than the total mass 
of the structure, however, the optically probed effective 
mass can have a strong dependence on the position and 
profile of the optical beam. We estimate the optically 
probed effective mass of the cantilever in our experiment 
using a combination of measurements and finite element 
analysis (FEA). Using the established values for the rele- 
vant elastic constants averaged over the crystalline multi- 
layer (Cii = 119.6, Ci2 = 55.5, C44 = 59.1 GPa) and the 
average material density 4476 kg/m^, the lateral geome- 
try of the FEA-simulated resonator is adjusted until min- 
imal error is found between the measured and simulated 
eigenfrequencies for the first four out-of-plane mechani- 
cal modes, see FigQa). (Note that the lowest frequency 
vibrational mode for our cantilever is an in-plane mode 
as the cantilever used is slightly thicker than wide.) A 
mean discrepancy between the measured and simulated 
frequencies of 6.1% was obtained by reducing the feature 
linewidth by 0.875 fim with respect to the lithographic 
mask. Note that the thickness of the free-standing mir- 
ror material was not used as a fitting parameter as it was 
accurately determined from the reflectance spectrum of 
the mirror [1 and found to be 6.88 fim. Once the geome- 
try is determined, the effective mass is calculated via the 
volume integral 



pfffdxdydz {u'^ 



(3) 



Here, p is the material density; and w are the dis- 

placements of the body along the x, and z directions, 
respectively, and the optically-probed displacement V is 
the overlap between the mechanical deflection and the 
optical Gaussian intensity profile, i.e. 



V 



dx dy w{x, y,z = 0) exp 



r 



2r2 



, (4) 



where ro is the standard deviation of the Gaussian op- 
tical intensity profile and the coordinate axis used for 
X, y and z has its origin in the center of the cantilever 
head, see Fig. ^h). The anti-reflection coated fiber fo- 
cuser used in our experiment provides an optical beam 
diameter (4ro) of 10.6 /im, which is much smaller than 
the nominal cantilever head diameter of 100 /im (as fabri- 
cated diameter of 98.25 /im). Thus, for the fundamental 
out-of-plane mode, there is only a weak optical beam po- 
sition and width dependence on the effective mass. (In 
this case, the effective mass is approximately equal to 



the intrinsic modal mass.) For the higher order modes of 
the structure, however, lateral displacement of the beam 
leads to a rapid change in the effective mass. In order 
to minimize the contribution from these higher mechani- 
cal modes it is necessary to carefully position the optical 
beam. Assuming careful alignment, the geometry of our 
mechanical structure is such that the contributions from 
higher order modes are further suppressed as the effec- 
tive mass rapidly increases with mode number since the 
majority of the deflection is due to the support beam for 
the fundamental mode rather than in the head for higher 
harmonics. Indeed, the unconditional RMS amplitudes 
of modes #4, #8, and #10 are 2.4%, 0.4%, and 0.1% 
that of mode #2, respectively. 





(941) 984 Hz 

#4: 3.68 

(11437) 10 975 Hz 
#8: 12.23 /ig^^^^^^^^^-^ 
^35 695) 33 215 Hz 
#10: 41.04 /ig 

(72 452) 66 863 Hz 



FIG. 4: (a) Finite element simulations of the vibrational 
modes (#2, 4, 8, 10) that have the four lowest optically probed 
effective masses. The displacements shown are exaggerated 
(blue to red indicates increasing displacement) and the outline 
indicates the cantilever rest position. The optically probed 
effective mass is given next to each mode number and under- 
neath are the simulated mechanical frequencies in brackets 
and the measured frequencies, (b) Schematic of the cantilever 
and the focussed signal beam. The coordinate axis used when 
determining the optically probed effective mass is in the cen- 
ter of the cantilever head on the mirror surface. 



Calibration Procedure 

We have used a two-step calibration procedure to de- 
termine the proportionality between the pulsed homo- 
dyne measurement outcomes and the mechanical dis- 
placement. During this procedure the signal beam is fo- 
cused onto the chip edge, i.e. a rigid unpatterned part 
of mirror material adjacent to the mechanical resonator, 
to prevent mechanical motion contributing to the sig- 
nal. First we calibrate the displacement of a piezoelec- 
tric actuator, which our fabricated structure containing 
the mechanical oscillator is placed upon, in response to a 
known drive voltage. We then drive the piezo and record 
the pulse measurement outcomes during the controlled 
actuation in order to calibrate the pulsed interferometer. 
Each step is detailed below in the next two subsections, 
respectively. 



9 




0.5 1 1.5 2 

time [ms] 



FIG. 5: Plot of the homodyne output (red, left axis) during 
piezo drive (blue, right axis) to calibrate the piezo response. 
The sinusoidal drive generates a peak-to-peak piezo scan of 
one half of an optical wavelength that was determined by 
adjusting the drive amplitude until the turning points in the 
homodyne output coincide (dashed line). 



was one half of the optical wavelength. This can be done 
precisely as the difference current output of the interfer- 
ometer has separate turning points occurring at the same 
level for this modulation depth and is then proportional 
to cos [00 + 7rsina;t], see Fig.jsj here 0o is the (unlocked) 
slowly varying phase in the interferometer and uj is the 
piezo drive angular frequency. As (pQ slowly changes this 
merely shifts the level of the turning points. In our exper- 
iment we used a drive frequency of 1.06 kHz and exploited 
a resonance of the piezo to achieve a peak-to-peak scan 
of 532 nm using 4.6 Vpp. 



Pulse Calibration 



Piezo Calibration 

To calibrate the piezoelectric actuator we applied a si- 
nusoidal drive voltage and used a continuous signal beam 
to monitor the piezo motion. The frequency of the drive 
was chosen such that the piezo mechanical response was 
either in or out of phase with the drive voltage. (Ex- 
perimentally, care was needed to find a suitable drive 
frequency as the piezo does not have a flat spectral re- 
sponse.) During this procedure, the phase between the 
signal and LO beams does not require locking and the 
piezo drive was at a higher frequency than the phase noise 
components in the interferometer. We then adjusted the 
drive amplitude such that the peak-to-peak piezo motion 



Using the same piezo drive frequency as above, and us- 
ing the piezo actuator calibration value (meters per Volt) 
obtained, the actuator was scanned with a reduced am- 
plitude so that the optical phase shifts are small. (It was 
verified that the piezo responds linearly with the applied 
Voltage over our range of interest.) Then, during the 
piezo scan, pulsed position measurements are performed 
and both the voltage applied to the piezo at the time 
of the measurement and the pulsed measurement out- 
comes are recorded. The proportionality between these 
recorded values is used to obtain the outcome per meter 
calibration. This calibration value is optical amplitude 
dependent and had to be measured for several optical 
amplitudes for the measurement shown in Fig. [2|f) of 
the main text. 
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